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Abstract 


An  explicit  expression  is  provided  for  the  characteristic  polynomial  of  a  matrix  M 
of  the  form 


M  =  D- 


0  ab'^  \ 
ba^  0  )  ’ 


(1) 


where  Z)  is  a  diagonal  matrix,  and  a  and  b  are  column  vectors.  Also,  an  explicit 
expression  is  provided  for  the  matrix  of  normalized  eigenvectors  of  M,  in  terms  of  the 
roots  of  the  characteristic  polynomial  (he.,  in  terms  of  the  eigenvalues  of  M). 


1  A  Lemma,  a  Remark,  and  an  Observation 

The  following  lemma  is  verified  by  substituting  into  the  left  hand  side  of  (7)  the  definitions 
of  P  in  (6)  and  U  in  (9)-(16),  and  simplifying  the  result  using  (4).  See  [2]  for  similar  results, 
and  [3]  and  [1]  for  applications. 

Lemma  1  Suppose  that  m  and  n  are  positive  integers,  a  =  (oq,  Oi,  . . . ,  0^-2,  and 

b  =  {bo,  bi,  . ,  bn-2,  bn-i)'^  are  real  vectors,  and  do,  di,  . ..,  dm+n-2,  dm+n-i  and  Xo,  Xi, 
. . . ,  Xm+n-2,  Xm+n-i  are  real  numbers  such  that 


Xj  7^  dk 

for  any  j,k  (j,k  =  0,  1,  . . . ,  m  +  n  —  2,  m  +  n  —  1), 

Xj  7^  Afc 


when  j  7^  k,  and 

{ak?  \  {b,?  \  ^  ^ 

\  fc=0  ~  )  \fc=0  —  Xj  J 

(with  j  =  0,  1,  ...,m  +  n  —  2,  m  +  n  —  1). 
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Suppose  further  that  D  is  the  diagonal  {m  +  n)  x  (m  +  n)  matrix  defined  by  the  formula 


(  do 

0 

0  \ 

0 

di 

dm+n—2  0 

V  0 

0  dra+n—1  ) 

and  P  is  the  (m  +  n)  x  (m  +  n)  matrix  defined  by  the  formula 

(  0  ab^\ 

^~[ba^  0  J  ’ 

where  0  denotes  matrices  consisting  entirely  of  zeroes. 

Then, 

{D-P)U  =  UK, 

where  A  is  the  diagonal  (m  +  n)  x  (m  +  n)  matrix  defined  by  the  formula 


(  Ao 

0 

0  \ 

0 

Ai 

Am+n— 2  0 

V  0 

0  ^m+n—l  ) 

and  U  is  the  orthogonal  {m  +  n)  x  {m  +  n)  matrix  defined  by  the  formula 

f  AVR  \ 


BWS  /  • 


In  (9),  A  is  the  diagonal  m  x  m  matrix  defined  by  the  formula 


(  ao 

0 

0  \ 

0 

ai 

O-m-2  0 

V  0 

0  1  / 

B  is  the  diagonal  n  x  n  matrix  defined  by  the  formula 


(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 
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V  is  the  m  X  {m  +  n)  matrix  with  entry  Vj^k  defined  by  the  formula 


Vj,k  = 


di  —  Xi 


(12) 


(with  j  =  0,  1,  ...,m  —  2,m  —  l;k  =  0,  1,  ...,m  +  n  —  2,m  +  n—l),W  is  the  nx  (m  +  n) 
matrix  with  entry  Wj^k  defined  by  the  formula 


Wxk  = 


—  At 


(13) 


(with  j  =  0,  1,  . . . ,  n  —  2,  n  —  1;  k  =  t),  1,  ...,m  +  n  —  2,  m  +  n  —  1),  S  is  the  diagonal 
(m+n)  X  (m+n)  matrix  with  the  diagonal  entries  So,o,  ■  ■  ■ ,  Sm+n-2,m+n-2,  Sm+n-i,m+n-i 
defined  by  the  formula 


\ 


m—1 

E 

k=0 


dk 


A.. 


n—1 

+  E 

fc=0 


^m-\-k 


A. 


(14) 


and  R  is  the  diagonal  {m  +  n)  x  (m  +  n)  matrix  with  the  diagonal  entries  Rofi,  Ri,i,  ■  ■  ■ , 
Rm+n—2,m+n—2j  Rm+n—i,m+n—i  defined  by  the  formula 


^j,j 


(15) 


In  (If)  and  (15),  Cq,  Ci,  . . . ,  Cm+n-2,  Cm+n-i  (ire  the  real  numbers  defined  by  the  formula 


n—1 

=  E 


{bkf 


fc=0  llm+k  -^j 


(16) 


Remark  2  The  equation  (4)  is  equivalent  to  the  characteristic  (secular)  equation 

det|Aj/-(T)-P)|  =  0  (17) 

for  the  eigenvalues  Xj  (with  j  =  0,  1,  . . . ,  m  +  n  —  2,  m  +  n  —  1)  of  the  matrix  D  —  P. 


Observation  3  The  upper  bloek  AVR  of  the  matrix  U  defined  in  (9)  has  the  form  of  a 
diagonal  matrix  (A)  times  a  matrix  of  inverse  differenees  (V )  times  another  diagonal  matrix 
(R).  The  lower  bloek  BW S  of  the  matrix  U  defined  in  (9)  also  has  the  form  of  a  diagonal 
matrix  (B)  times  a  matrix  of  inverse  differenees  (W )  times  another  diagonal  matrix  (S). 
Therefore,  there  exists  an  algorithm  which  applies  such  an  N  x  N  matrix  U  (or  its  adjoint) 
to  an  arbitrary  real  vector  of  length  N  in  O  {N  log(l/e))  operations,  where  e  is  the  precision 
of  computations  (see  [3]). 


2  Acknowledgements 

The  author  would  like  to  thank  V.  Rokhlin  for  many  useful  discussions  and  proofreading. 


3 


References 


[1]  S.  Chandrasekaran  and  M.  Gu,  A  divide- and- conquer  algorithm  for  the  eigende- 
composition  of  symmetric  block- diagonal  plus  semiseparahle  matrices,  Numerische  Math- 
ematik,  96  (2004),  pp.  723-731. 

[2]  M.  Gu  AND  S.  G.  Eisenstat,  A  stable  and  efficient  algorithm  for  the  rank-1  modifica¬ 
tion  of  the  symmetric  eiqenproblem,  SIAM  Journal  of  Matrix  Analysis  and  Applications, 
15  (1994),  pp.  1266-1276. 

[3]  - ,  A  divide- and- conquer  algorithm  for  the  symmetric  tridiagonal  eigenproblem,  SIAM 

Journal  on  Matrix  Analysis  and  Applications,  16  (1995),  pp.  172-191. 


4 


